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Abstract 

Pulsed  power  gas  switches  operate  under  extreme 
conditions  holding  off  voltages  at  millions  of  volts  and 
greater,  conducting  currents  at  100’s  of  kiloamps,  and  can 
reliably  trigger  with  less  than  10ns  of  jitter  at  optimal 
operating  conditions.  However  as  Martin  and  Braginskii 
had  noted  the  impedance  of  these  switches  is  dynamic  and 
a  function  of  many  parameters  such  as  gas  pressure, 
current,  and  the  number  of  channels.  In  systems  such  as 
Protogen  or  Ursa  Minor  many  gas  switches  are  utilized  to 
enable  optimal  performance  and  pulse  shaping.  A 
challenge  with  either  system  is  that,  if  the  switches  are 
tightly  coupled  and  one  switch  triggers  early,  then  the 
fields  on  the  second  switch  change  dramatically  and  can 
move  the  second  switch  into  a  new  operating  regime  that 
will  result  in  a  different  dynamic  resistance.  This  paper 
presents  a  least  squares  approach  to  extracting  the  Martin 
and  Braginskii  switch  resistance  with  field  dependent 
carrier  recombination  from  measured  data. 

I.  INTRODUCTION 

Switches  are  a  critical  component  in  new  pulsed  power 
technologies  that  are  seeking  to  enhance  system 
performance  and  attain  extreme  changes  in  operational 
regimes.  Furthermore,  new  technologies  utilizing  more 
modular  approaches  such  as  Protogen  [1]  or  linear 
transformer  driver  (LTD)  [2]  based  systems  (Ursa  Minor 
[3],  Mykonos  [4],  or  the  next  generation  of  the  Z  machine 
[5])  are  heavily  dependent  upon  the  synchronized 
operation  of  hundreds  to  thousands  of  gas  insulated 
pulsed  power  switches. 

Simultaneous  operation  of  many  switches  is  a 
challenging  problem  in  itself  with  sensitivity  to  variations 
in  the  manufacturing  process,  trigger  resister  impedance 
variations  due  to  manufacturing  and  aging,  as  well  as 
variations  in  trigger  generator  performance  when  multiple 
trigger  generators  are  operated  in  parallel.  These 
variations  result  in  changes  to  differences  in  electric  field 
stresses  from  switch  to  switch  which  can  dramatically 
impact  jitter  and  runtime.  This  impacts  efficiency  and 
energy  density  of  pulsed  power  systems. 

Output  pulse  shaping  through  precision  control  of 
switch  closing  time  is  a  technique  that  continues  to  be 
refined  for  dynamic  materials  research  [1].  Systems  such 
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Figure  1.  The  switch  circuit  and  idealized  components. 
All  are  assumed  constant  except  Rd. 


as  Protogen  [1]  and  the  Z  machine  [6],  [7]  have 
demonstrated  such  capabilities  however  changes  in  field 
level  for  switches  that  are  delayed  in  time  exacerbates  the 
increase  of  jitter  and  run  time.  Techniques  that  avoid 
using  switch  control  to  define  output  pulse  shape  have 
been  demonstrated  in  systems  such  as  Veloce  [8]  but  with 
limited  success  in  obtaining  a  wide  range  of  pulse  shapes 
with  variations  not  only  in  rise  time  but  also  in  the  shape 
of  the  peak  and  the  falling  edge  of  the  waveforms. 

This  paper  presents  results  on  characterization  of  a 
modified  High  Current  Electronics  Institute  (HCEI)  LTD 
200  kV  gas  switch  [9]  over  a  large  range  of  operating 
conditions.  The  equivalent  circuit  diagram  representing 
the  test  setup  for  collecting  the  switch  data  is  provided  in 
Figure  1.  A  photograph  of  the  two  stage  Marx  consisting 
of  the  switch,  two  100  kV,  80  nF  capacitors,  and  a  low 
impedance  load  (CVR)  is  in  Figure  2.  In  this  fixture  both 
the  voltage  across  the  capacitor  and  the  current  through 
the  switch  are  measured. 


Figure  2.  The  brick  configuration  of  the  test  circuit. 
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In  previous  papers,  data  from  this  circuit  was  fit  using 
the  Braginskii  [10]  -  Martin  [11]  model  with  (1)  genetic 
optimization  [12],  and  (2)  piecewise  constant  resistance 
least  squares  fits  (LSFs)  to  a  damped  harmonic  oscillator 
solution  [13].  The  genetic  optimization  gave  reasonable 
fits,  but  was  time  intensive  (hours  on  one  PC).  Piecewise 
LSFs  were  much  faster  (~30s).  With  this  approach,  the 
time  period  for  the  current  was  divided  into  quarter  period 
intervals  (—150  ns),  and  each  interval  was  fit  using 
constant  switch  resistance  and  inductance.  The  best  fit 
switch  resistance  was  then  plotted  for  each  interval, 
displaying  a  slowly  varying  resistance.  However,  the 
piecewise  fit  provided  very  limited  time  resolution  of  the 
dynamic  switch  resistance,  because  the  intervals  must  be 
relatively  large  to  get  reasonable  fits  to  the  current.  Since 
switch  resistance  changes  very  rapidly  as  the  switch 
closes  (at  early  time),  a  higher  resolution  solution  was 
desired.  Both  approaches  could  also  be  improved  by 
including  more  realistic  plasma  behavior  such  as:  field- 
dependent  carrier  generation,  carrier  recombination,  and 
carrier  density  non-uniformities.  In  this  paper,  a  least 
squares  fit  is  made  to  the  current  in  the  differential 
equation  of  state  (DES)  using  a  modified  Braginskii- 
Martin  model  for  the  switch  resistance  that  includes 
electric  field  dependent  carrier  generation  and  carrier 
recombination  in  a  cylindrical  switch  plasma  channel  with 
changing  diameter.  Then  a  numerical  solution  is 
calculated  for  the  current  using  the  fit  coefficients  for  the 
DES.  This  provides  a  single,  high  resolution,  time- 
dependent  solution  to  the  current  with  dynamic  switch 
resistance  that  decreases  rapidly  at  early  times  and 
increases  slowly  at  late  times  as  the  current  decays. 

II.  SWITCH  MODEL 


The  switch  model  chosen  to  represent  the  performance 
of  the  HCEI  switch  is  based  on  spark  channel  work  by 
Braginskii  [10]  and  switch  loss  work  by  Martin  [11].  It 
uses  a  time  varying  resistance  to  describe  the  energy  loss 
in  the  switch.  This  particular  model  has  the  advantage  of 
being  simplistic  in  form  thereby  enhancing  the  efficiency 
of  system  current  shaping  optimization  algorithms.  In  this 
model  the  radius  of  the  spark  channel,  a,  is  a  function  of 
the  current  flowing  through  the  channel  as  defined  by 
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where  /  is  the  magnitude  of  the  current  in  the  channel,  p0 
is  the  gas  density,  a  is  the  conductivity  of  the  gas,  and 
is  set  to  a  constant  value  of  4.5  as  defined  in  [11].  The 
resistance  of  a  single  channel 

Rcn=-^~2  (2) 
a  n  az 

is  then  based  on  the  radius,  where  lch  is  the  channel 
length.  From  Rch  the  total  HCEI  dynamic  switch 
resistance  is  calculated  as 
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where  nch  is  the  total  number  of  spark  channels  formed. 

If  switch  inductance  and  all  other  circuit/component 
values  are  considered  constant  except  the  dynamic  switch 
resistance,  Rd ,  then  the  DES  describing  the  voltage  across 
the  ideal  switch  represented  in  Figure  1  is 


^0  -  (I'c  +  +  ( Rc  +  Rs  +  RdV>  t))I  +  c 
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where  subscripts  C  and  S  stand  for  fixed  circuit/ 
component  and  switch  values,  V0  is  the  initial  voltage 
across  the  capacitor,  and  t0  is  the  switching  start  time. 

The  dynamic  resistance  is  approximated  by 
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where  n0  and  n(t)  are  the  initial  and  temporal  ion  number 
densities  with  <r(t)  =  n(t)eii  and  the  constant 

(6) 


r  (n0eM)2/3(47r/(PoD)1/3 
The  Braginskii-Martin  models  assume  temporally 
constant  conductivity  in  the  gas.  In  this  paper,  the  ion 
density  in  the  switch  channel  is  modeled  as  a  function  of 
the  electric  field  governed  by  the  rate  equation  with  terms 
for  avalanche  generation  and  carrier  recombination  with 
the  following  constraints. 
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E  (t)  is  the  magnitude  of  the  electric  field  in  the  channel 
E0  is  the  field  at  which  h=0  (balanced) 

Tr  is  the  zero  field  exponential  recombination  time 
Ta  is  the  high  field  exponential  carrier  generation  time 
Ea  is  the  field  at  which  /  (Ea)  =  ^/ja 

The  general  solution  for  the  ion  number  density  is 

n(t)/n0  =  eSto  f{E^dt'  (8) 

The  electric  field  dependence  of  the  rate  equation, 
was  modeled  four  ways:  linear,  quadratic, 
exponential,  and  hyperbolic  (Figure  3).  Only  the 
quadratic  and  hyperbolic  models  could  fit  estimated 
values  for  the  three  constraints  in  equation  (7)  (Tr~10‘6s, 
ra~10'9s,  £o~2.1x105V,  £a~2.5xl05V).  To  do  this,  the 
quadratic  model  went  to  a  negative  minimum  between 
E  =  0  and  E0  which  is  unphysical.  Hence  the  most 
appropriate  model  was  the  hyperbolic. 
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Figure  3.  Rate  equation  models  for  a  200  kV  switch. 

The  electric  field  in  the  switch  channel  is  approximated 
with  the  state  equation  4  using  an  initial  estimate  for  the 
fixed  switch  inductance  and  resistance. 

-  Vo -(Lc+Ls0)i-(Rc+Rs0)l-S  Idt'/C  „„ 
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The  DES  was  fit  with  the  initial  estimates.  Then  the 
results  were  substituted  and  the  fit  was  iterated.  Stable 
values  were  reached  in  2  iterations. 

III.  ANALYSIS  OF  RESULTS 

The  current  in  the  switch  circuit  (Figure  4)  was 
measured  with  a  high  bandwidth  current  viewing  resistor 
(CVR).  Digital  noise  in  these  measurements  was  reduced 
by  filtering  the  data  with  a  low  pass  filter  (Butterworth 
30-50  MHz  passband-stopband).  Substituting  equations 
11,  9,  8,  and  5  into  the  DES  (4),  produces  a  system  of 
approximations  for  each  current  data  point  with  three 
unknown  linear  coefficients:  Ls,  Rs,  ar.  The  sum  of  the 
squared  residuals  in  these  approximations  due  to  noise 
can  be  minimized  by  solving  for  unique  values  of  these 
coefficients  using  least  squares  fitting  (LSF). 

The  fixed  component  values,  C,  Rc ,  and  Lc,  were 
measured  using  a  precision  impedance  analyzer.  The 
values  used  for  these  calculations  were  37.52  nF,  200  mQ, 
and  200  nH.  Other  constants  in  these  equations,  nch,  xr, 
Ta,  Ea,  Eq,  t0,  tf0,  Ibo/  were  estimated  and  adjusted 
manually  to  minimize  the  residuals  in  the  fit.  (t0  is  the 
switching  start  time,  tf0  is  the  filter  start  time,  Ib0  is  the 
current  baseline  offset.)  Including  any  of  these  constants 
as  fitting  parameters  will  require  a  non-linear  FSF,  due  to 
the  dynamic  resistance  dependence  on  current  (equation 
5).  The  convergence  of  such  non-linear  FSF  may  be 
tested  in  future  analyses. 

A  200  kV  switching  example  is  shown  in  Figure  4.  The 
raw  (red  dots)  and  filtered  data  (red  line)  are  compared  to 
solutions  to  the  DES  for  3  different  FSFs:  (1)  fixed 
resistance,  (2)  constant  resistivity,  and  (3)  dynamic  radius 
and  resistivity.  Forcing  a  fixed  switch  resistance  (green), 


Figure  4.  The  circuit  current  for  a  200  kV  example  is 
measured  with  a  high  bandwidth  CVR  and  shown  full 
scale  (top),  early  time  enlarged  (middle),  and  late  time 
enlarged  (bottom)  as  raw  data  (red  dots),  filtered  data 
(red),  and  the  solutions  to  the  DES  with:  (1)  fixed  switch 
resistance  (green),  (2)  dynamic  switch  resistance,  constant 
ion  number  density  (black),  and  (3)  dynamic  switch 
resistance,  hyperbolic  ion  number  density  (blue). 
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in  the  fit,  results  in  a  current  that  starts  abruptly  from  zero 
with  a  magnitude  that  is  too  large  at  the  beginning  and  the 
end  of  the  pulse.  Allowing  a  dynamic  switch  resistance 
and  constant  ion  number  density  (black),  in  the  fit, 
improves  agreement  near  switching  initiation,  but  has  too 
large  a  magnitude  near  the  end  of  the  pulse.  The  best  fit 
was  achieved  with  dynamic  resistance  and  the  hyperbolic 
ion  number  density  (blue).  Both  fits  using  dynamic 
resistance  start  smoothly  from  zero  current. 

The  fit  coefficients  and  the  dynamic  resistance  from  the 
example  shown  in  Figure  4  are  shown  in  Figure  5  with 
their  95%  confidence  intervals  of  uncertainty.  The 
enlargement  (bottom)  shows  the  current  periodicity  in  the 
dynamic  resistance.  Note  the  initial  fall  in  resistance  the 
model  predicts  with  the  sudden  increase  in  current.  The 
analytic  models  for  channel  diameter  and  the  ion  number 


Figure  5.  The  dynamic  resistance  was  calculated  using 
the  modified  Braginskii-Martin  model  with  hyperbolic 
field-dependent  ion  number  density  to  fit  the  DES. 


density  provide  solutions  with  smooth  current  turn-on  as 
the  switch  closes.  This  solution  was  insensitive  to  the 
avalanche  carrier  generation  time,  as  the  effective 
bandwidth,  -500  MHz,  and  oscillation  frequency,  ~1.7 
MHz,  of  the  data  are  too  small  to  measure  sub¬ 
nanosecond  phenomena  associated  with  sudden  switch 


Time  (us) 


Figure  6.  The  normalized  ion  number  density  calculated 
for  the  example  shown  in  the  previous  figures.  While  the 
fit  to  the  data  was  very  sensitive  to  the  recombination 
time,  Tr,  it  was  insensitive  to  the  avalanche  carrier 
generation  time,  ra. 

closure.  The  late  time  resistance  increase  is  very  sensitive 
to  carrier  recombination  as  the  field  in  the  switch  channel 
decreases.  The  normalized  ion  number  density  for  this 
example  is  shown  in  Figure  6. 

Indications  for  the  quality  of  the  fit  to  the  DES  voltage 
and  the  quality  of  the  DES  solution  for  current  are 
displayed  in  Figure  7.  The  voltage  residuals  settle  quickly 
from  the  initial  switch  closure  to  less  than  2%  of  the 
initial  voltage,  200  kV.  Residuals  for  the  current  using 
dynamic  resistance  with  hyperbolic  ion  number  density 
are  shown  relative  to  the  raw  data  (red  dots)  and  the 
filtered  data  (red  line).  The  filtered  residuals  also  settle 
quickly  to  less  than  2%  of  the  initial  current  (70  kA).  At 
late  time,  they  are  small  compared  to  the  noise  in  the  raw 
data.  Also  plotted  are  the  residuals  for  the  current  using 
fixed  resistance  relative  to  the  filtered  data  (green  line). 
These  match  the  dynamic  resistance  results  in  the  middle 
of  the  time  range,  but  are  much  larger  near  the  beginning 
and  end  of  switching.  The  adjusted  R-square  (normalized 
measure  of  the  quality  of  fit)  for  this  solution  based  on  the 
current  residuals  was  0.995  (a  perfect  fit  would  be  1). 

To  assess  the  reproducibility  of  these  results,  6  shots  at 
200  kV  initial  switch  voltage  and  the  same  pressure  were 
analyzed.  The  switch  resistance,  Rd  +  Rs,  for  each  shot  is 
compared  in  Figure  8.  The  example  in  the  preceding 
discussion  and  figures  is  the  red  curve  (shot  43)  in  these 
plots.  The  coefficients  in  the  state  equation,  inductance, 
Ls,  and  dynamic  resistance,  ar,  of  the  switches  are 
provided  in  table  I  along  with  the  total  resistance  minima, 
Rm,  and  start  times,  t0.  The  top  plot  in  this  figure  displays 
the  dynamic  resistance  over  the  full  time  range,  -8  ps. 

The  other  plots  in  Figure  8  are  enlargements  which 
show  an  early  time  interval  (middle)  and  an  interval 
containing  the  resistance  minima  (bottom).  The  initial 
range  of  variation  in  time  (0.4  ns)  is  due  to  scope  shot-to- 
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Figure  7.  Residuals  in  voltage  from  the  DES  fit  and  in 
current  from  the  DES  solutions  are  plotted  for  the 
example  in  the  previous  figures. 

shot  sampling  synchronization,  as  the  sampling  period  is 
1.3  ns.  Timing  jitter  was  removed  by  basing  start  times 
(see  t0  in  Table  I)  on  the  maximum  of  the  2nd  derivative  of 
the  filtered  data.  The  linear  segments  in  this  plot  are  due 
to  the  linear  extrapolation  of  the  current  derivative  back  in 
time  to  the  time  when  the  smoothed  current  was  zero. 
The  standard  deviation  in  the  resistance  minima  (2mQ)  is 
very  small  when  compared  to  the  overall  range  in  the 
switch  resistance  calculation,  10  kQ  (top),  and  the 
amplitude  noise  in  this  data,  ~  2  kA  out  of  1 60  kA  peak- 
to-peak  (Figure  4). 

The  variation  in  the  total  switch  resistance  as  a  function 
of  the  initial  switch  voltage  at  the  pressure  to  hold-off  200 
kV  is  shown  in  Figure  9  and  Table  II.  The  minimum  total 
switch  resistance,  Rm ,  (bottom)  and  run  time  (delay  from 
trigger,  t0  in  Table  II)  decrease  with  increasing  charge 
voltage.  Initial  range  of  time  and  amplitude  variation  is 
similar  to  the  variation  shown  Figure  8  (middle). 


Time  (ns) 


-8  -7  -6  -5  -4  -3 

Time  (ns) 


Figure  8.  Comparison  of  the  total  switch  resistance,  Rs  + 
Rd ,  calculated  for  6  shots  at  200  kV.  The  top  plot  shows  the 
full  scale  of  calculation.  The  middle  plot  is  an  enlargement 
of  the  onset  of  switching.  The  bottom  is  an  enlargement  of 
the  interval  containing  the  resistance  minima. 
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Table  I.  The  fit  coefficients  for  six  shots  at  200  kV. 


shot 

43 

44 

45 

46 

47 

48 

mean 

std.  dev. 

Rm(mQ) 

46 

44.3 

49.8 

48 

47.7 

48.7 

47.4 

1.96 

Ls(nH) 

26.1 

23.8 

28.5 

21.4 

25.9 

28.7 

25.7 

2.80 

ar(xl06) 

7.46 

7.78 

7.54 

7.71 

6.32 

6.93 

7.29 

0.560 

to  (ns) 

884 

886 

885 

886 

887 

886 

886 

1.13 

Table  II.  The  fit  coefficients  for  a  range  of  initial  charge  voltages. 


Voltage  (kV) 

100 

120 

140 

160 

180 

200 

mean 

std.  dev. 

Rm(mQ) 

137 

107 

69.2 

67.2 

50.8 

46 

79.5 

35.4 

Ls(nH) 

55.8 

48.4 

37.6 

31.4 

23.5 

26.1 

37.1 

12.8 

2s  N 

O 

X 

a 

7.52 

7.24 

10.2 

7.13 

6.38 

7.46 

7.66 

1.32 

t0  (ns) 

989 

943 

911 

902 

894 

884 

921 

39.3 

Figure  9.  Total  switch  resistance  versus  charge  voltage. 


IV.  CONCLUSIONS 

The  Braginskii  and  Martin  switch  resistance  with 
field-dependent  carrier  recombination  was  extracted  from 
measured  data  by  utilizing  least  squares  fitting  and 
solving  the  differential  equation  of  state.  Modeling  the 


switch  performance  as  a  function  of  current  magnitude, 
spark  channel  radius,  avalanche  generation,  and  carrier 
recombination  improved  the  prediction  of  switch 
performance  from  the  closing  through  the  opening 
process.  Both  the  dynamic  switch  resistance  and  run  time 
were  observed  to  decrease  with  increases  in  the  initial 
charge  voltage. 
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